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ABSTRACT 

Aims. Detailed information about the magnetic geometry, atmospheric abundances and radial velocity variations has been obtained for 
the magnetic standard star HD 94660 based on high-dispersion spectroscopic and spectropolarimetric observations from the UVES, 
HARPSpol and ESPaDOnS instruments. 

Methods. We perform a detailed chemical abundance analysis using the spectrum synthesis code zeeman for a total of 17 elements. 
Using both line-of-sight and surface magnetic field measurements, we derive a simple magnetic field model that consists of dipole, 
quadrupole and octupole components. 

Results. The observed magnetic field variations of HD 94660 are complex and suggest an inhomogeneous distribution of chemical 
elements over the stellar surface. This inhomogeneity is not reflected in the abundance analysis, from which all available spectra 
are modelled, but only a mean abundance is reported for each element. The derived abundances are mostly non-solar, with striking 
overabundances of Fe-peak and rare-earth elements. Of note are the clear signatures of vertical chemical stratification throughout the 
stellar atmosphere, most notably for the Fe-peak elements. We also report on the detection of radial velocity variations with a total 
range of 35 km s“* in the spectra of HD 94660. A preliminary analysis shows the most likely period of these variations to be of order 
840 d and, based on the derived orbital parameters of this star, suggests the first detection of a massive compact companion for a main 
sequence magnetic star. 

Conclusions. HD 94660 exhibits interestingly complex magnetic field variations and remarkable radial velocity variations. Long term 
monitoring is necessary to provide further constraints on the nature of these radial velocity variations. Detection of a companion will 
help establish the role of binarity in the origin of magnetism in stars with radiative envelopes. 
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1. Introduction 

Of order 10% of main sequence A and B stars are magnetic. For 
the majority of these stars, their magnetic properties are studied 
using circular polarisation spectra, which provide a measure of 
the mean magnetic field projected along the line-of-sight 
A fraction of these stars rotate slowly and/or have strong enough 
magnetic fields that Zeeman splitting is obs erved in indiv idual 
lines of the intensity (Stokes I) spectra (see iBailevL I20l4 and 
references therein). The magnitude of this splitting provides a 
measure of the mean magnetic field modulus at the stellar sur¬ 
face ((B)). These magnetic A and B stars exhibit anomalous at¬ 
mospheric abundances and are referred to as the chemically pe¬ 
culiar A- type (Ap) stars. These stars often have (Bj.) in excess 


* Based in part on our own observations made with the European 
Southern Observatory (ESO) telescopes under the ESO programme 
093.D-0367(A) and programmes 076.D-0169(A), 088.D-0066(A), 
087.D-0771(A), 084.D-0338(A), 083.D-1000(A) and 60.A-9036(A), 
obtained from the ESO/ST-ECF Science Archive Facility. It is also 
based in part on observations carried out at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
of Canada, the Institut National des Science de I’Univers of the Centre 
National de la Recherche Scientifique of France and the University of 
Hawaii. 


of about 1 kG and ( B) of order several thousand s to tens of thou¬ 
sands of gauss (e.g. lDonati & Landstreelll2009ll . 

The chemical peculiarities of Ap stars can be quite striking. 
For example, the Fe-peak element Cr may be as much as 10^ 
times overabundant compared to the Sun. Still more impressive 
are the abundances of rare-earth elements, which are commonly 
in excess of solar values by as much as 10^ times. 

Ap stars may be highly variable, with their magnetic field 
strengths and spectral line strengths and shapes varying with 
the rotation period of the star. T his variability is best explained 
via the rigid rotator model (see IStibbsL Il950l) . In this model, 
the line-of-sight and magnetic axis are at angles i and /3 to 
the rotation axis, respectively. Therefore, different magnetic 
field measurements throughout the rotation cycle of the star are 
the result of observing the field at different orientations and 
since chemical elements are distributed non-uniformly and non- 
axisymmetrica lly over the stellar su rface, spectrum variability is 
also observed (IRvabchikoval ll 99 ll) . 

HD 94660 ( = HR 4263) is a bright (V = 6.11) chemically 
peculiar magnetic Ap star that is commonly used as a magnetic 
standard to test polarimetric syste ms in the southern hemisph ere. 
The field was first discovered bv iBorra & Landstreell (1 197.51) us¬ 
ing an Ha magnetograph; they reported a value for the line-of- 
sight magnetic field of (Bj) = -3300 + 510 G. HD 94660 is a 
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Table 1. Summary of the stellar and magnetic properties of 
HD 94660. 


Spectral type 

AOp EuSiCr 

Renson et al., 1991 

Teff (K) 

11300 ±400 

This paper 

log g (cgs) 

4.18 ± 0.20 

This paper 

R(Ro) 

2.53 ± 0.37 

This paper 

vsini (kms“‘) 

<2 

This paper 

P(d) 

2800 ± 250 

Landstreet et al.. 2014 

log(L*/Le) 

2.02 ±0.10 

This paper 

M(Mo) 

3.0 ± 0.20 

This paper 

n (milliarcsec) 

6.67 ± 0.80 

van Leeuwen 2007 

Bd (G) 

-7500 

This paper 

Bq (G) 

-2000 

This paper 

B„c, (G) 

7500 

This paper 

in 

16 

This paper 

pn 

30 

This paper 


sharp-lined star with clearly res o lved Zeeman splitting in sev¬ 
eral spectra lines (lMathvslll990h . lBohlender et al.l(ll993h report 
a projected rotational velocity of v sin i < 6 km s“’ and a roughly 
constant line-of-sight magnetic field strength of (B^) - -2520 G, 
based on measurements using YiB- A period of ro tation of order 
2700 d w as first proposed by iHensbergel (Il993h and later dis¬ 
cussed b v iMathv s et al. ( 1997h with respect to (B) data. More 


m 

recently, ILandstreet et al.l (120141) studied the field variations of 
HD 94660 from 17 FORSl observations taken over a 6 year pe¬ 
riod from 2002 to 2008. They find a peak-to-peak variation in 
(Bj) of about 800 G, ranging from about -2700 and -1900 G. 
From these variations, a rotation period of 2800 + 250 d is de- 
du ced, which agrees with previo us determinations. 

ILandstreet & Mathvsl (|2000|) were the first to report clear 
(Bj) variations (between ab out -1800 and - 2 100 G ) which, 
together with (B) data from iMathvs & Hubrid (Il997h . enabled 
them to model the magnetic field of this star using a colin- 
ear multipole expansion. They adopted a model that consists of 
dipole, quadrupole and octupole components with surface po¬ 
lar field strengths of -8400, 2700 and 6900 G, respectively. As 
seen above, the value of (B^) is always negative, indicating that 
i+/3< 90° i.e. only the negative magnetic pole is observed. This 
is confirmed by this model where the determined values for i and 
y6 are 5° and 47°. We point out that, in general, this model does 
provide a rough first approximation to the field variations. 

In this paper, we discuss efforts to model the magnetic field 
and chemical abundances of many elements based on high- 
dispersion, polarimetric spectra. The following section discusses 
the observations. Sect. 3 outlines the derived physical param¬ 
eters; Sect. 4 discusses the magnetic field measurements and 
model; Sect. 5 and 6 describe the modelling technique and abun¬ 
dance analysis, respectively; Sect. 7 reports detected radial ve¬ 
locity variations; and Sect. 8 summarises the results of the paper. 


2. Observations 

We have acquired one spectropolarimetric observation of 
HD 94660 and have retrieved an additional seven archival spec¬ 
tra that were utilised in this study, all of which were taken 
with HARPSpol. This is a high-resolution (B ^ 115 000), 
cross-dispersed echelle spectropolarimeter that covers a spectral 
range between 3780 - 6910 A and is mounted on the European 
Southern Observatory (ESO) 3.6-m telescope located at La Silla. 
The high signal-to-noise ratio (SNR) observations consist of 6 
circularly polarised Stokes V spectra (with the 01-04-2012 ob¬ 


servation also including linear polarisation in both Stokes Q and 
U), and two unpolarised Stokes I spectra, all acquired between 
May 2009 and April 2014. Each polarimetric observation (ex¬ 
cept for the observation taken on 31-05-2009) was obtained by 
acquiring four successive individual spectra with the quarter- 
wave plate rotate d in such a way to acqui re the desired Stokes 
spectrum (see e.g. lRusomarov et al.Ll2013L for further details). 

The HARPSpol spectra were reduced using a modified 
version of the r e duce package (iPiskunov & Valeiitil 120021 : 
iMakaganiuk et all 1201 ih . Wavelength calibration was per¬ 
formed using the spectrum of a ThAr calibration lamp and 
then corrected to the heliocentric rest frame. Normalisation of 
the spectra was achieved by first dividing the spectra by the 
optimally-extracted spectrum of the flat-held to correct for the 
blaze shape and fringing. The resulting spectra were then cor¬ 
rected by the response function derived from observations of 
the Sun. The last step involved fitting a smooth, slowly-varying 
function to the envelope of the entire spectrum. The final output 
is a set of continuum-normalised Stokes I spectra, the Stokes 
parameter of interest (V, Q, U) and a diagnostic null spectrum 
that is calculated in such a way that the polarisation cancels out, 
which often allows us to identify spurious signals that are present 
in the processed data. The spectrum that was acquired on 31- 
05-2009 only completed half of the full polarimetric sequence, 
which enables the cancellation of first order (linear) wavelength 
drift, but does not correct for second order (quadratic) drift 
and does not allow the computation of a diagnostic null spec¬ 
trum. We also note that the retarder angles listed in the obser¬ 
vations obtained on 05-01-2010 were inconsistent with the usu¬ 
ally adopted values for obtaining Stokes V measurements and 
are likely erroneous. We therefore proceeded to reduce the data 
assuming the usual retarder angles for the sequence of obser¬ 
vations, and then verified that the resulting polarised and unpo¬ 
larised spectra were in good agreement with the other observa¬ 
tions. 

In addition to the HARPSpol observations, we identified an¬ 
other archival, circularly polarised, high-resolution (B ^ 65 000) 
spectropolarimetric observation acquired with ESPaDOnS on 
January 9, 2006. ESPaDOnS, which is mounted at the 3.6 m 
Canada-Erance-Hawaii Telescope (CEHT), is also a bench- 
mounted, cross-dispersed echelle spectropolarimeter, which 
covers a broader spectral range compared to HARPSpol of 3690 
- 10 481 A. The polarimetric observation was obtained by tak¬ 
ing a sequence of four sub-exposures with different positions of 
the Eresnel Rhomb to acquire a singl e Stokes V spectrum , ac¬ 
cording to the procedure described bv iDonati et alJ (Il997l) . The 
ESPaDOnS spectra were processed using the automated reduc- 
tion package libre-e sprit, following the double-ratio procedure 
dPoriati et al.Lll99’^ . 

Lastly, we also found three nights of archival UVES data 
consisting of a total of 11 spectra. UVES is a cross-dispersed 
spectrograph mounted on the ESO 8.2-m Very Large Telescope 
(VLT) located at Paranal. It has both a blue and red arm offering 
different spectral resolutions and wavelength coverage. The blue 
arm offers B up to about 80000 and a spectral range of 3100 - 
4900 A, whereas the red arm covers 4800 - 10200 A with B 
up to about 110000. Tablel^summarises our entire collection of 
spectra for HD 94660. 
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Table 2. Log of available spectra for HD 94660. Successive columns list the instrument, date and JD of observation, exposure time, 
estimated signal-to-noise ratio (SNR) per 1.8 km s“^ velocity bin at 5000 A, spectral resolution (R) and wavelength coverage. 


Instrument 

Date 

(DD-MM-YYYY) 

JD 

(2450000+) 

^exp (s) 

SNR 

R 

A ( A) 

UVES 

01-05-2001 

2031.464 

145 

254 

65 030 

3043-3916 



2032.464 

70 

252 

74 450 

4726 - 6808 



2031.466 

70 

266 

74 450 

4726 - 6808 

UVES 

01-08-2001 

2038.441 

100 

379 

65 030 

3731 -4999 



2038.443 

100 

383 

65 030 

3731 -4999 



2038.441 

100 

101 

74 450 

6650 - 10 426 



2038.443 

100 

101 

74 450 

6650 - 10 426 

UVES 

03-12-2005 

3707.846 

200 

233 

71 050 

3044-3917 



3707.841 

200 

316 

71 050 

3281 -4563 



3707.846 

200 

675 

107 200 

4726 - 6835 



3707.841 

200 

283 

107 200 

5708 - 9464 

ESPaDOnS 

09-01-2006 

3745.167 

1200 

270 

65 000 

3690- 10 481 

HARPSpol 

24-05-2009 

4975.546 

600 

311 

115 000 

3780-6910 

HARPSpol 

25-05-2009 

4976.536 

1200 

398 

115 000 

3780-6910 

HARPSpol 

31-05-2009 

4982.605 

1200 

476 

115 000 

3780-6910 

HARPSpol 

05-01-2010 

5201.833 

808 

397 

115 000 

3780-6910 

HARPSpol 

19-05-2011 

5701.450 

800 

392 

115 000 

3780-6910 

HARPSpol 

20-05-2011 

5702.449 

800 

427 

115 000 

3780-6910 

HARPSpol 

01-04-2012 

6018.559 

1000 

670 

115 000 

3780-6910 

HARPSpol 

28-04-2014 

6775.610 

360 

320 

115 000 

3780-6910 


3. Physical Parameters 

3.1. Effective Temperature and Gravity 

Geneva and Stromgren uvbyP photometry are available for 
HD 94660 and both were utilised to determine the effective tem¬ 
perature Leff and gravity log g for the star. For the Geneva pho¬ 
tometr y, we used the fortran c ode developed by iKunzli et alJ 
(Il997h . A modified version of the lNaniwotzki et al] ( 1993h cod^ 
that corre cts the effective t emper ature to the Ap temperature 
scale (see iLandstreet et all l2007l for a complete discussion), 
was used for the Stromgren photometry. 

From the Geneva and Stromgren photometric systems, we 
found Teff = 11500 K, log g = 4.10 and - 11100 K, log g - 
4.26 respectively. We adopt the mean of these two sets of values 
for our analysis with Teff = 11300 + 400 K and log g - 4.18 
+ 0.2, with the uncertainties estimated from the scatter between 
the measurements, and by taking into account the intrinsic un¬ 
certainties in computing these parameters for Ap stars. These 
parameters are used for our abundance analysis (see Sect. 6). 

3.2. Luminosity, Steiiar Radius and Mass 

Ivan LeeuwenI (l2007h reports a Hipparcos parallax of 6.67 + 0.80 
milliarcseconds. From this value, a distance to HD 94660 of 
about 150 pc is deduced. With a well determined distance, an ap¬ 
propriate bolometric correction for Ap stars can be use d to deter¬ 
mine the stellar luminosity (see lLandstreet et ^l2007ll . which is 
found to be log L/Lq = 1.97+0.25 with uncertainties propagated 
in the usual way. 

Based on the and luminosity determinations, it is 
straightforward to compute a stellar radius of R = 2.53 ±0.37 Rq. 
By further comparing the position of HD 94660 to theoreti¬ 
cal evolutionary tracks (iGirardi et al.L 1200^ in an HR diagram, 
we are able to estimate the star’s evolutionary mass. Using the 
adopted uncertainties in T^ff and L, we proceed by comparing the 
star’s position to multiple evolutionary tracks of varying masses. 
In this manner, we estimate that HD 94660 has a mass of about 


3.0 + 0.20 Mq. This is in agreement to the mass of 3.51 + 
0.64 Mq that we can estimate from the stellar radius and photo¬ 
metrically determined log g. Further, the position in the HR dia¬ 
gram suggests that HD 94660 has completed le ss than half of its 
main s equence lifetime. This is consistent with ILandstreet et"^ 
(l2008ll . who find strong magnetic fields only in young Ap stars. 

4. Magnetic Field 


4.1. Longitudinal Magnetic Fieid Measurements 

For each of the Stokes V spectra we measured the mean, surface- 
averaged longitudinal magnetic fie ld ((B,,)) from line- averaged 
Least-Squares Deconvolved (LSD; IPonati et all 1 19971) line pro- 
hles. The technique involves obtaining mean line profiles by 
combining all spectral lines in a given line list (normally metallic 
and He lines). The result is a much higher SNR with the ability 
to detect weaker Zeeman signatures due to magnetic helds. 
wa s measured using the first-order moment meth od discussed 
by Rees & Semel J 19791). and as implemented by iDonati et al.l 
(11997 ) and Wade et al.l J 2 OOOI) according to the equation: 


-2.14x10" [(v - vo)F(v)dv 

(Bz) = -^- —r -• 1 ) 

4zc /[l-/(v)]dv 

In this equation, v is the velocity within the LSD prohle, V{v) 
is the continuum-normalised Stokes I prohle and I{v) is the 
continuum-normalised intensity prohle. The wavelength A (in 
nm) and Lande factor z correspond to the weighting factors used 
in the calculation of the LSD prohles (500-nm and 1.2, respec¬ 
tively). 

The LSD prohle s were extracted using the iLSD code 
of iKochukhov et al.l (I 2 OIOI) . As input the code requires a 
line mask that w as extracted from the Vienna Atomic Line 
Database (VALD; iKunka et^ l2000t iRvabchikova et al.L Il997t 
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Fig. 1. (Bj> and (B) variations of HD 94660. The solid black lines 
are the adopted magnetic field model variations from Sect. 4.3. 
Top: (Bj.) LSD measurements using the full mask (black circles) 
and Nd mask (green triangles). Note that the open and filled 
symbols denote measurements from the ESPaDOnS spectrum 
and HARPSpol spectra, resp ectively. Balmer line measurements 
from iLandstreet et~an (|2014|) are the filled red squares. Bottom: 
Phased (B) measurements from Fe ii 6149 (blue squares), Cr ii 
6833 (red triangles) and Nd iii 6145 (green circles). The dotted 
green line is the best-fit sinusoidal variations to the Nd measure¬ 
ments. 



0 2 4 6 8 10 

Normalised EQW 


Fig. 2. The variations in magnetic field strength (Bj) versus 
EQW for individual spectral lines. Different symbols denote dif¬ 
ferent elements, while different colours are used to represent 
different Lande factor ranges. Note that the absolute value of 
the field strength is shown and each measurement has been nor¬ 
malised as discussed in the text. The dashed line is a linear fit to 
the data. 


of interest and the other containing all other lines. Both masks 
are used simultaneously as input and allow us to extract a repre¬ 
sentative mean, unblended profile of the element of interest. All 
LSD profiles were extracted onto a 1.8 km s ' velocity grid and 
uncertainties were computed by propagating the uncertainties in 
the final LSD profiles. 

Our final measurements are listed in Table |3] The top panel 
of Figure [Uplots the (Bj) measurements against rotational phase 
for the full metallic and Nd masks as well as the magnetic field 
model derive d in Sect. 4.3. Also show n are the Balmer line mea¬ 
surements of iLandstreet et al] (l2014l) . Although the HARPSpol 
and ESPaDOnS magnetic data have not been intercalibrated for 
this star, we note that the consistency between the ESPaDOnS 
and HARPSpol measurements is satisfactory. Furthermore, pre¬ 
vious studies have shown a good agreement between the po- 
larimetric spectra acq uired with ESPaDOnS and HARPSpol 
(iPiskunov et al.Ll201 Ih . 


iPiskunov et al.L Il995t iKupka et al.L 1 19991) for the spectral range 
covered by HARPS, using the mean abundances determined in 
this work and discussed in Sect. |6l This linelist was used for 
both the HARPS and ESPaDOnS spectra to provide a consis¬ 
tent (Bj) measurement between the two datasets. We then pro¬ 
ceeded to remove all lines that were blended with lines not used 
in this analysis, such as broad hydrogen lines or lines blended 
with strong telluric absorption bands. As a final step we then 
automatically adjust the line depths from their theoretical pre¬ 
dictions to provide a best fit to the observed Stokes I spectrum, 
while also removing lines that are poorly fit (such as those lines 
that show very strong Zeeman splitting). This final mask was 
then used to extract mean line profiles for all spectra that we 
label as ‘full mask’. We also extracted mean line profiles of indi¬ 
vidual elements for Si, Ti, Cr, Fe, and Nd, using the multi-profile 
capabilities of iLSD. In each case we provided two input masks 
based on our final full mask, one made entirely of the element 


4.1.1. Measurements from individual spectral lines 


The extremely sharp lines and strong magnetic field result in 
clean, strong Stokes V signatures for individual spectral lines. 
This allows us to test the multi-line technique of LSD (see 
Sect. 4.1) by measuring (Bj.) from individual lines of elements. 
We performed these measurements for a total of about 80 un¬ 
blended lines which included the elements Si, Ti, Cr, Fe and Nd 
(the same elements for which we performed LSD to compute the 
field from individual elements). In general, the measurements 
agree with the values we derive and present in Table [3] but with 
large scatter between different lines of an individual element. 
In an attempt to understand this scatter, we plot the measured 
(Bj) value against the equivalent width (EQW) for all the spec¬ 
tral lines (see Fig. EJ. To account for any variation in the field 
measurements due to the different phases, all values have been 
normalised by the (B^) and EQW value for the Fe 1 4982 line oh- 









































J. D. Bailey et al.: The Magnetic Standard HD 94660 


5 


tained for a given night. Note that the absolute value of the field 
strength is shown for clarity. 

The results show a clear trend, namely a stronger field 
strength for lines with smaller EQW and a weaker field strength 
for lines with larger EQW. We note that almost all of the outliers 
from this trend were lines with very small Lande factors. These 
results are difficult to interpret. On one hand, this result could re¬ 
flect differential desaturation i.e. that strong lines should be more 
sensitive to the desaturation of the line due to magnetic splitting 
and so the relative change in the EQW of that line should be 
greater than for a weak line in the presence of a strong magnetic 
field. We attempted to test this hypothesis by producing a syn¬ 
thetic spectrum using zeeman (see Sect. 5), which includes all 
of the measured spectral lines. Because the field strength is kept 
constant for a given phase in zeeman for all lines, any systematic 
change in the measured {B^) value as a function of the EQW can 
be attributed to differential saturation. Unfortunately, the results 
we derived from this test are also difficult to interpret, but it ap¬ 
pears that the general trend of the synthetic tests do not agree 
with the results from the observations. Another possible expla¬ 
nation caused by differential saturation is that this result could 
simply reflect small systematic changes in the measured centre- 
of-gravity of Stokes V (the numerator in Eqn.[T]i. This would pre¬ 
dominantly affect strong lines with complex splitting patterns, 
resulting in a blended profile that does not have the same centre- 
of-gravity as the unsaturated line with the same splitting pattern. 

On the other hand, this result could be unrelated to saturation 
effects and could be a measure of the change of the decreasing 
field strength with increasing vertical height in the atmosphere: 
the cores of weak lines are formed deeper in the atmosphere and 
should have intrinsically stronger fields relative to strong lines 
with cores formed higher. Using our zeeman model to estimate 
the depth of the atmosphere (~0.6% of the stellar radius), we 
compute an expected change of the order of ~4% between the 
measured field strength of lines formed at the top and bottom of 
the atmosphere. Eig.|2] shows a typical variation in the strength 
of the magnetic field of order 30% from weak to strong lines 
(closer to 50% if you consider the extreme values). Therefore, it 
seems unlikely that this effect is important and the nature of this 
trend is still unclear. 


4.2. Surface Magnetic Field Measurements 


The sharp-lined features and high resolving powers of the 
HARPSpol, UVES and ESPaDOnS instruments allowed us to 
make measurements of the mean field moduli ((B)) from ob¬ 
served Zeeman splitting in the lines of Nd iii /16145, Ee ii 46149 
and Cr ii 46833. Notable exceptions are for the one ESPaDOnS 
spectrum and two nights of UVES spectra. The lower resolving 
power of ESPaDOnS compared to HARPSpol did not allow us to 
see sufficient splitting in Nd iii, whereas insufficient wavelength 
coverage for the UVES spectra did not allow for measurements 
of splitting in Cr ii (on 01-05-2011) or Nd iii and Ee ii (on 01-08- 
2011 ). 

To obtain values of (B), the atomic data from the VALD 
database was used and the field strength calculated from the 
equation 

“ 4.67 X 10-13^2^’ 


where A4 denotes the shift of the cr components from the zero 
field wavelength, 4o is the rest wavelength in A and z is the ef¬ 
fective Lande factor of the line (1.00, 1.34 and 1.33 for 46145, 
6149 and 6833, respectively). The separation of the elements of 
the Zeeman split lines that determine A4 was measured by fitting 
Gaussians to each component using the splat function in IRAE. 
The uncertainties were estimated by considering the dispersion 
between a series of meas urements of ( B) fro m multiple Gaussian 
fits to each component. iBailev et’aP (1201 ih note that measure¬ 
ments of (B) from Eqn. |2] are very sensitive to small changes 
in A4. Because of the high resolution and essentially zero rota¬ 
tion rate of HD 94660, the dispersion in our measurements of 
A4 is less than about 0.005 A in most cases. Table|4]summarises 
these (B) measurements and the bottom panel of Eig.[T]show the 
phased rotational variations of the surface fields from the three 
elements. 

The scatter of (B) values around the mean curve confirms 
that our uncertainty estimates are reasonable. The consistency of 
measurements from each individual element taken on successive 
nights (24-05-2009 and 25-05-2009) and the agreement, within 
estimated uncertainties, between (B) measurements of the Ee- 
peak elements Cr ii and Ee ii (except for a single measurement) 
further suggests that the quoted uncertainties are realistic. The 
surface field appears to vary sinusoidally; however, the average 
surface field from the Fe-peak elements is of order 400 G larger 
than that of Nd iii. Further, the peak-to-peak amplitude of the 
variations is twice as large for Nd iii (of order 800 G) compared 
to the Fe-peak elements (of order 400 G). The fact that different 
elements sample the field in significantly different ways suggest 
that the abundance distributions of different elements, in addition 
to the field structure, are inhomogeneous over the stellar surface. 
This is further supported by the discrepant (B^) measurements of 
different elements discussed above and listed in Table [3] 

4.3. Magnetic Field Model 

Given the long rotation period of HD 94660, we are fortunate 
to have sufficient phase coverage for both (B^) and (B) in order 
to derive a simple magnetic field model. In general, the FORS1 
data show that the Balmer line and metallic line (B^) measure- 
ments agree closely, ex cept for one discrepant metal (B^) value 
(iLandstreet et al.Lr2014l) . However, a large discrepancy exists be¬ 
tween all FORS 1 measurements and those from high-resolution 
spectra. In particular, the FORS 1 data show a much larger range 
in (Bj). This can be seen in the top panel of Fig. [1] where the 
FORS 1 Balmer line measurements vary by order 800 G whereas 
the full LSD mask from the high resolution spectra vary by about 
300 G. W e opt for using the Balm er line (B^) variations as dis¬ 
cussed by ILandstreet et all (l2014t) to derive a field geometry. We 
point out that this is the main reason for the discrepancy be- 
tween our field struc t ure (s ee below) and the one derived by 
ILandstreet & Mathvsl(l200(^ . 

_ The FORTRAN program fldsrch.f (ILandstreet & Mathvsl 

I 2 OOOI) derives a magnetic field geometry based on the observed 
(Bj) and (B) variations. It produces a simple co-axial multipole 
expansion that consists of dipole (Bd), quadrupole (Bq) and oc- 
tupole (Boa) components with the angles between the line-of- 
sight and rotation axis (i) and the magnetic field and rotation axes 

o 

(j8), specified. For HD 94660, a new field geometry with / = 16 , 
^ = 30°, Bd = -7500 G, Bq = -2000 G, and Bo^t = 7500 G 
is found. The phased (B^) and (B) variations using our adopted 
geometry are shown in Fig. [T] and are in good agreement with 
the longitudinal field variations measured by ILandstreet et al.l 
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Table 3. Log of the line-of-sight surface field, {B.), measurements of HD 94660. Listed are the instrument, date and JD of observa¬ 
tions, phase and {B^) measurements from the full stellar mask and only lines of Si, Ti, Cr, Fe and Nd. 


Instrument 

Date 

(DD-MM-YYYY) 

JD 

(2450000-t-) 

Phase 

Full Mask 

Si 

<Sz> (G) 

Ti Cr 

Fe 

Nd 

ESPaDOnS 

09-01-2006 

3745.167 

0.051 

-2223 ± 17 

-2343 ± 77 

-1963 ±61 

-2174 ±25 

-2217 ±19 

-2601 ± 70 

HARPSpol 

31-05-2009 

4982.605 

0.493 

-2287 ±14 

-2231 ±71 

-2341 ±51 

-2114 ±22 

-2329 ±16 

-2381 ±64 

HARPSpol 

05-01-2010 

5201.833 

0.571 

-2413 ±14 

-2460 ± 70 

-2444 ± 50 

-2273 ± 23 

-2480 ±17 

-2628 ± 58 

HARPSpol 

19-05-2011 

5701.450 

0.750 

-2304 ±14 

-2515 ±68 

-2281 ±51 

-2106 ±22 

-2348 ±17 

-2739 ± 53 

HARPSpol 

20-05-2011 

5702.449 

0.750 

-2303 ±14 

-2502 ± 68 

-2276 ±51 

-2092 ± 22 

-2351±17 

-2726 ± 52 

HARPSpol 

01-04-2012 

6018.559 

0.863 

-2252 ±13 

-2399 ± 65 

-2352 ± 48 

-2151 ±22 

-2280 ±16 

-2625 ± 49 

HARPSpol 

28-04-2014 

6775.610 

0.133 

-2151±14 

-2305 ±71 

-1643 ±48 

-2149 ±21 

-2200 ±17 

-2372 ± 61 


Table 4. Log of surface field, {B), measurements of HD 94660. For each spectra, the instrument, date and JD of observation, phase 
and {B) measurements from Nd iii 46145, Fe it 46149 and Cr it 46833 are listed. 


Instrument 

Date 

(DD-MM-YYYY) 

JD 

(2450000±) 

Phase 

NdiiiT6145 

<S> (G) 

Fe II T6149 

Cr II T6833 

UVES 

01-05-2001 

2031.464 

0.439 

6124 ± 184 

6292 ± 94 


UVES 

01-08-2001 

2038.441 

0.441 

_* 

_* 

6253 ± 120 

UVES 

03-12-2005 

3707.841 

0.038 

5954 ± 179 

6711 ± 105 

6253 ± 80 

ESPaDOnS 

09-01-2006 

3745.167 

0.051 

- 

6460 ± 84 

6315 ± 102 

HARPSpol 

24-05-2009 

4975.546 

0.490 

5877 ± 129 

6123 ± 60 

6208 ± 117 

HARPSpol 

25-05-2009 

4976.536 

0.490 

5868 ± 141 

6176 ± 63 

6191 ±77 

HARPSpol 

31-05-2009 

4982.605 

0.493 

5900 ± 113 

6122 ± 98 

6221 ± 119 

HARPSpol 

05-01-2010 

5201.833 

0.571 

5881 ± 114 

6093 ± 84 

6194 ± 147 

HARPSpol 

19-05-2011 

5701.450 

0.750 

5395 ±71 

6085 ± 64 

6044 ± 89 

HARPSpol 

20-05-2011 

5702.449 

0.750 

5549 ± 93 

6039 ± 60 

5993 ± 88 

HARPSpol 

01-04-2012 

6018.559 

0.863 

5708 ± 83 

6220 ±51 

6244 ± 78 

HARPSpol 

28-04-2014 

6775.610 

0.133 

6209 ± 85 

6418 ± 43 

6426 ± 60 


Notes. (*) UVES spectra do not include these spectral lines (see Table[2Jl 


( 12014 fromFORSl, as well as the measured Fe-peak surface 
field variations. 

Although our magnetic held model is based on the ht to 
the large amplitude {B^) FORSl data, it is nevertheless useful 
to com pare our geometry to the one from iLandstreet & Mathvsl 
(I 2 OOOI) because it is not entirely clear which, if either, is bet¬ 
ter to use. In general, our adopted held geometry is not dras¬ 
tically different than the one proposed by iLandstreet & Mathvsl 
(I 2 OOOI) (see Sect. 1). The angles i + (5 are roughly the same, al¬ 
though the individual values of these angles differ (recall that i 
and p can be interchanged). The values of B^ and Boct are sim¬ 
ilar; however, the relatively small Bq, although comparable in 
magnitude, is n egative in our derived geom etry and positive in 
the geometry of iLandstreet & Mathvsl (120001) . The circular polar¬ 
isation signatures in Stokes V are also well produced with both 
geometries. Fig. [3] compares the observed and synthetic Stokes 
I and V prohles, using both models, for Nd iii 46145 (note that 
only one phase is shown but the splitting is equally well repro¬ 
duced at all phases). It is obvious from this figure that both helds 
do an adequate job at reproducing the observed magnetic prop¬ 
erties in the spectrum of HD 94660. Based on Figures [3l it is 
unclear which model is best and both models are too coarse to 
describe the complex field variations over the visible stellar sur¬ 
face. Although our field model is preferred, we perform an abun¬ 
dance analysis using both geometries (see Sect. 5). 

ILandstreet & Mathvsl (l2000t) suggest that fi is likely smaller 
than i for a very large fraction of slow rotators, arguing that for 
long period stars (B > 25 d) it is equally likely for the rotation 
pole to be positioned anywhere on the visible hemisphere. For 
this reason, smaller values of i are statistically unlikely, since 


it would require i to be within a small region of the visible 
hemisphere. However, their geometries are based only on inten¬ 
sity and circular polarisation spectra which makes it difhcult to 
distinguish i and fi. Fortunately, we have available two obser¬ 
vations of the linear polarisation (Stokes Q and U) spectra of 
HD 94660 taken with the HARPSpol instrument on 01-04-2012 
(see Sect. 2). The linear polarisation data will allow us to dis¬ 
tinguish which of the two angles is larger. Figure |4] compares 
the Stokes Q signature for our adopted magnetic field model, 
assuming that i > /3 (top panel) and i < p (bottom panel). For 
i > jS, the model clearly predicts a Stokes Q signature that is 
too large compared to the observations. Alternatively, for i 
the model predicts a much more modest signature, in relative 
agreement to what is observed, although we admit that the model 
does not accurately represent the true variations in linear polar¬ 
isation for this star (we n ote that similar result s are a lso seen 
with Stokes U and with the lLandstreet & Mathvsl(l2000l) model). 
This hgure strongly suggests that, contrary to the assumption of 
ILandstreet & Mathvsl ( 20001) . for HD 94660 p is the larger of 
these two angles for both models. 


5. Spectrum Synthesis Program 

5.1. ZEEMAN.F 

The FORTRAN program zeeman.f is a spectrum synthesis program 
developed bv ILandstreet! (Il988h . zeeman.f performs line forma- 
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Fig. 3. A comparison of magnetic field models. In black are 
the observed Stokes I (top) and V (bottom) profiles for Nd iii 
/16145 for the 28-08-2014 HARPSpol spectrum. We overplot 
the computed line profiles for our adopted magnetic field ge- 
ometry (red dashed line) a nd the magnetic field geometry of 
iLandstreet & Mathvsl (l2000l) (blue dotted line). 



Fig. 4. A comparison of the Stokes Q signatures of synthetic 
models to the HARPSpol observation on 01-04-2012. In black 
is the observed Stokes Q spectrum. Overplotted is our adopted 
magnetic field model (red dashed line). The top panel assumes 
i > p, whereas the bottom assumes i < p. 


tion and radiative transfer in a magnetic field. The input mag¬ 
netic field geometry consists of the axisymmetric superposition 
of a dipole, quadrupole and octupole with specified values of the 
angles i and p. zeeman.f interpolates an appropriate atmospheric 
model from a grid of solar abundance ATLAS9 models based 
upon the Tetr and log g supplied. The atomic data needed for line 
synthesis are taken from the VALD database, zeeman.f assumes 
a homogeneous abundance distribution vertically throughout the 
stellar atmosphere and computes all four Stokes parameters, zee¬ 
man.f allows up to six abundance rings that are symmetric about 
the magnetic axis to be specified on the stellar surface, with each 
ring having equal extent in colatitude. Within each ring, a uni- 


Table 5. The abundances are measured with respect to H, and are 
tabulated with their associated uncertainties. For reference, the 
solar abundance ratio and the number of spectral lines modelled 
are also shown. 


Element 

logfAx/Nn) 

Solar 

# Lines 

He 

< -2.20 

-1.07 

2 

0 

-3.77 ± 0.20 

-3.31 

3 

Mg 

-5.03 ±0.15 

-4.40 

1 

Si 

-3.21 ±0.30 

-4.49 

5 

Ca 

-6.19 ±0.33 

-5.66 

1 

Ti 

-5.55 ±0.50 

-7.05 

7 

Cr 

-3.86 ±0.30 

-6.36 

30 

Mn 

-5.67 ± 0.20 

-6.57 

2 

Fe 

-3.10 ±0.50 

-4.50 

28 

Co 

-4.51 ±0.25 

-7.01 

8 

Ni 

-5.24 ± 0.20 

-5.78 

1 

Sr 

-7.78 ± 0.40 

-9.13 

2 

La 

-7.68 ±0.16 

-10.9 

2 

Ce 

-6.86 ± 0.20 

-10.42 

3 

Pr 

-6.80 ± 0.25 

-11.28 

4 

Nd 

-6.54 ± 0.30 

-10.58 

9 

Eu 

-7.49 ± 0.20 

-11.48 

3 


form abundance is assumed. For each modelled spectrum, a best- 
fit V sin i and radial velocity vr are calculated. 

5.2. Choice of magnetic fieid modei 

One of the major difficulties in performing the abundance anal¬ 
ysis of a magnetic star is the inclusion of the effects of the mag¬ 
netic field on spectral line formation. This is, in part, due to the 
limited number of tools that include this physics, but also be¬ 
cause of the inherent ambiguity in selecting an appropriate mag¬ 
netic field model. Section 4 highlighted this discrepancy with 
two apparently adequate field models to explain the observed 
Stokes I and V profiles of HD 94660. Therefore, we compared 
the abundances derived from these two magnetic field geome¬ 
tries. We also tested the abundances with a simple dipolar model 
in which the value of was chosen to be approximately three 
times the largest observed {Bp value (this leads to roughly the 
same extrema in the observed {Bp and (B) values). 

We note that the agreement between the different models is 
impressive. For any given set of fundamental parameters (i.e. Tefr 
and log g), the difference in the derived abundances between the 
three models are generally less than about 0.1 dex and no worse 
than about 0.2 dex. It is encouraging that the choice of model 
does not drastically influence the final abundances. In fact, a 
very coarse model, consisting only of a dipole, can give accurate 
abundances, which suggests that the final abundances may not be 
very sensitive to the difference between these simple models and 
the real field distribution. Apparently, the inclusion a magnetic 
field model that roughly explains the {Bp and (B) variations is 
sufficient to provide an acc urate analysis. This lends support to 
studies such as the one by iBailev et al.l (1201 4l) for which sim¬ 
ple dipolar models are used to quantify the evolution of atmo¬ 
spheric abundances in Ap stars. Since the abundances depend 
little on the magnetic field model, the following section reports 
only on abundances found using the field geometry we present 
in Sect. 4.3. 
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Fig. 5. Example of three synthesised spectral windows for the 
28-04-2014 HARPSpol spectrum of HD 94660. The observed 
spectrum is in black (solid line) and the model is in red (dashed 
line). The bottom panel illustrates the vertical stratification of Fe 
in the atmosphere of HD 94660, where strong lines of Fe it and 
weak lines of Fe i are fit well with the abundance of Table|5] but 
the weaker lines of Fe it require an enhanced abundance to be 
well modelled. 


6. Elemental Abundances 

Insignificant variability in the strengths and shapes of most spec¬ 
tral lines is observed among the twelve spectra available to 
model for HD 94660. For lines that exhibit strong Zeeman sig¬ 
natures (for example Nd iii at 6145 A) differences can be seen in 
the line structure (i.e. relative strengths of the cr and tt compo¬ 
nents) and, in some cases, marginal changes in the line strength. 
However, line depth changes are limited to less than about 0.05 
in intensity. Therefore, although all spectra are modelled, we re¬ 
port only a mean abundance for each element. Uncertainties are 
determined in two ways, depending upon the number of spectral 
lines available to model. For elements with multiple lines, the 
uncertainties are estimated from the observed scatter between 
the best-fit abundances. When only one (or few) lines are avail¬ 
able to model, we estimate uncertainties by changing the best-fit 
abundance for each element until the reduced deviates from 
best-fit models by about 1 (i.e.;^^^ = Xhest uncertainty is 

then the difference between the two abundances and corresponds 
approximately to a 1 cr uncertainty. Table|5]lists the abundances 
of each element, as well as the total number of lines modelled 
for each e lement. For reference , the solar abundance ratios are 
provided (lAsnlund et al.L [200^ . Fig. |5] provides an example of 
the quality of fit for three spectral windows for the HARPSpol 
observation taken on April 28, 2014. Below, each element is dis¬ 
cussed individually. 


6.1. Helium 

Many lines of helium exist to derive an abundance. He i /14471 
and 45876 are predicted to be the strongest helium lines in the 
spectrum, but they are not detected in the available spectra. An 
upper limit for helium is derived from these regions that is at 
least a factor of 12 below the solar abundance. As expected of a 
magnetic Ap star, HD 94660 is He-weak. 

6.2. Oxygen 

The O I multiplet at 6155-56-58 A is present in all spectra and is 
used to derive a mean abundance. The ESPaDOnS spectrum also 
includes the O i lines at 7771-74-75 A, however, this triplet suf¬ 
fers from non-FTE effects and is therefore not considered. The 
adopted abundance is about 2.5 times less than in the Sun. 

6.3. Magnesium 

The only line suitable for modelling is Mg ii at 4481 A. At all 
phases, this line is reasonably well fit with an abundance that is 
a factor of 4 below the solar ratio. 


6.4. Siiicon 

Several clean lines of Si ii exist throughout the spectrum includ¬ 
ing strong lines at 5041 and 5055-56 A, as well as the weaker 
line at 4621 A. When fit simultaneously, 45041, 5055-56 re¬ 
quire an abundance that is about 0.5 dex less than that of 44621. 
Table |5]reports the average of these abundances. 

We also found a line of Si iii at 4552 A. This line requires an 
abundance that is more than 10 times that of the mean abundance 
from the Si ii lines. This is exactl y what was reported for a sam - 
ple of mid to late B-type stars bv iBailev & Fandstreell (l2013al) . 
who argue that the discrepancy between the abundances derived 
from the first and second ionisation states of silicon is likely the 
result of strong vertical stratification in the stellar atmosphere. 
The different abundance values found using stronger and weaker 
lines of Si ii, and the line of Si iii, indicate that stratification is 
also likely in HD 94660. 

6.5. Caicium 

Possible useful lines of calcium at 8498, 8542 and 8662 A are 
present in the ESPaDOnS spectrum, however, these are blended 
with the Paschen lines, which cannot be calculated correctly with 
Zeeman at present. Therefore, only Ca ii at 3933 A was used to 
derive the final abundance of calcium. We note that distinctly 
different abundances are required to adequately fit the wings and 
core of this line (see iBabell 1 19921 where this effect is first ex¬ 
plained for Ca). An abundance derived from the wings of this 
line, which are formed deeper in the atmosphere, suggest an 
abundance that is approximately solar. On the other hand, the 
core (formed higher up in the atmosphere) requires an abun¬ 
dance that is of order 0.8 dex below the solar abundance. This is 
a symptom of vertical stratification and it would appear calcium 
is strongly stratified throughout the atmosphere of HD 94660. 
The mean abundance of calcium is reported in Table |5] 
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6.6. Titanium 

For deriving the mean abundance of titanium, there are several 
lines to choose throughout the blue spectrum. We have modelled 
the lines of Ti ii at 4533,4549,4563, 4568 and 4571 A, and also 
at 4798 and 4805 A. Each set of lines is synthesised separately 
and within each set the lines are modelled simultaneously. The 
final abundance is the mean of these two values (see Table |5]l 
The abundance that fits well the weaker lines of 44568 and 
4798, as well as the weak lines that are blended with other Fe- 
peak elements at 44533 and 4549, is systematically too large for 
the stronger lines at 445 63, 4571 and 4805. T his is a symptom 
of vertical stratification (iBagnulo et al.L ItOOH) . Nevertheless, Ti 
is clearly overabundant compared to the Sun by about a factor of 
30. 

6.7. Chromium 

Unlike the other Fe-peak elements, the abundance derived for Cr 
is not as drastically dependent upon the lines modelled, which 
suggests that Cr may not be significantly stratified throughout 
the stellar atmosphere. This is illustrated in Fig. |5] where most 
lines of Cr ii (both strong and weak), as well as two lines of 
Cr I, are fit well with the same uniform abundance. Note that 
HD 94660 is relatively hot for the presence of neutral atoms, 
however, these lines are unambiguously detected in all modelled 
spectra and are predicted to be visible by zeeman.f. Furthermore, 
tests of the depth of core formation for these lines indicate that 
they are formed in the upper part of the stellar atmosphere, where 
the temperature is necessarily cooler and consistent with the 
presence of neutral chromium. However, over 30 lines of Cr ii 
were tested throughout the spectrum with this uniform abun¬ 
dance, and a notable fraction of weaker lines require a larger 
abundance than is recorded in Table |5] to be adequately mod¬ 
elled. Although apparently not as significant as other elements, 
stratification may still be present for Cr. The final abundance of 
Cr was found by simultaneously fitting Cr ii lines at 4531, 4539, 
4555, 4558, 4565, 4587,4588, 4590 and 4592 A. At all phases, 
Cr is modelled well with an abundance that is of order 300 times 
larger than the solar ratio. 

6.8. Manganese 

The abundance of Mn was derived from the pair of Mn ii lines at 
6122 A. The lines are well modelled at all phases with a uniform 
abundance that is of order 8 times larger than the solar abun¬ 
dance. 

6.9. Iron 

One mean abundance is not adequate to satisfactorily model all 
the spectral lines of Fe. It is clear from Fig. |5] that the adopted 
abundance fits the weaker lines well, but is systematically too 
large for the stronger Fe ii lines. This is evident throughout the 
spectrum, with well over 20 Fe lines modelled with the uniform 
abundance listed in Table|5] This is symptomatic of vertical strat¬ 
ification in the atmosphere. Multiple lines of Fe ii are used to 
compute the mean abundance of iron and the larger estimated 
uncertainty is indicative of the discrepancies found when mod¬ 
elling each line. The weaker Fe ii lines at 5029, 5031, 5032, 
5034, 5035 and 5036 were fit simultaneously, and required an 
abundance larger by 0.6-0.7 dex than that deduced using strong 
lines such as 4583 and 5018 A. Three lines of Fe i at 5014, 6136 


and 6137 A were also modelled and required a similar abun¬ 
dance to what is required to fit the stronger Fe ii lines. It is clear 
that Fe is at least of order 30 times more abundant than in the 
Sun. 

6.10. Cobait 

The spectrum appears extraordinarily rich in Co ii. Several lines 
of Co II are scattered throughout the spectrum. A total of eight 
lines were modelled; 5016, 5017, 5023,5025, 5027,5129, 5135, 
5131 A. All lines of Co ii are reasonably well modelled with a 
uniform abundance and suggest an abundance that is at least 300 
times greater than that of the Sun. 

The detection of cobalt in the spectrum of a hot Ap star is rare 
(we note that i t is more common in the cooler roAp stars such 
as 10 Aquilae: iNesvacil et alil20T3h . We are only aware of two 
other Co-strong hot magnetic Ap stars (HR 1094 an d HR 5049; 
iNielsen & Wahlgreiil 120^001: iDworetskv et al.L Il980h . The abun¬ 
dances for both these Co-strong stars were derived without in¬ 
cluding the effects of the magnetic field and suggest abundances 
of this element greater than 1000 times that of the Sun. It is un¬ 
clear how prevalent Co is in the atmospheres of magnetic Ap 
stars, which emphasises the importance of further detailed anal¬ 
yses of other magnetic Ap stars. 

6.11. Nicket 

Ni II at 4067 A was used to derive the mean abundance. The line 
is reasonably well modelled in all the spectra and Ni is appar¬ 
ently overabundant compared to the solar abundance ratio by a 
factor of about 3. 

6.12. Strontium 

Sr II lines at 4077 and 4215 A are present in the blue spec¬ 
trum at all phases. The accurate modelling of these lines depends 
strongly on the adopted abundance of Cr ii with which the two 
Sr lines are blended. Since the blending Cr lines are weak lines 
in the wings of the Sr lines, we use a somewhat larger abun¬ 
dance for Cr than the value listed in Table |5] when modelling the 
lines of strontium. In this way, we get more concordant results 
between these two lines of Sr ii, but the overall agreement is still 
poor. Nevertheless, it is clear that Sr is overabundant, probably 
by a factor of about 25 compared to the Sun. 

6.13. Lanthanum 

Two lines of La ii that are suitable for modelling were found in 
the spectrum at 4605 and 6126 A. The abundance was derived 
from the former and tested using the latter. This rare-earth ele¬ 
ment is over 2000 times more abundant than in the Sun. 

6.14. Cerium 

Several lines of Ce ii are available to model throughout the spec¬ 
trum including 4560,4562 and 4628 A. The final abundance was 
found from modelling 44628, which satisfactorily fits the other 
lines of Ce. The models suggest that this rare-earth element is 

more than 3000 times the solar abundance ratio. _ 

Inte restingly, unlike what was found bv iBailev & Landstreetl 
(1201 3bh for HD 147010, there is no discrepancy between the 
abundances derived from 4560-62 A compared to 4628 A and 
suggests that perhaps the discrepancy they report between these 
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lines for that star may not be due to inaccurate gf values, but 
instead possibly due to an unrecognised blend or blends in this 
more rapidly rotating magnetic star (which has a v sin i of about 
15 kms^'). 

6.15. Praseodymium 

Lines of Pr iii at 4625, 6160 and 6161 A are available to model 
in all spectra. For the ESPaDOnS spectrum, we were also able 
to model Pr iii at 7781 A. Little to no variation is observed in the 
strength of these spectral lines and a uniform abundance models 
the observed spectrum well at all phases. Similarly to the other 
rare-earth elements, Pr is dramatically more abundant than in the 
Sun, by a factor of order 3000. 

6.16. Neodymium 

Nd has the richest spectrum of all the rare-earth elements with 
multiple lines of Nd iii available for modelling; 4570, 4625, 
4627, 4911, 4912, 4914, 5050, 5127 and 6145 A. In general, 
abundances derived from each line agree well with one another 
and are reasonably well modelled at all phases. Nd has the high¬ 
est abundance of any rare-earth element studied here, being of 
order 10"^ times the solar value. 

6.17. Europium 

Two suitable lines of Eu ii are present for modelling; 4129 and 
6645 A. One line of Eu iii is also available at 6666 A, however, 
this line is badly blended with Fe-peak elements. This rare-earth 
element is about 1000 times more abundant than in the Sun. 

7. Radial Velocity Variations 

In the process of performing our spectroscopic analysis, it be¬ 
came evident that HD 94660 exhibits significant radial velocity 
variations. We thus report substantial radial velocity (RV) vari¬ 
ations in the magnetic standard star HD 94660, with a range in 
ve locities of order 35 k m s '. These variations were first reported 
bv iMathvs et al.l dl997h . who noted that the orbital period should 
not be more than about 2 years, significantly sh orter than the 
rotation period. Some years later, iMathvsl (l2013li determined a 
value of the orbital period for HD 94660 of 848.96 days. Our 
search for the best-fit period of these radial velocit y variations 
was c arried out using the Lomb-Scargle method dPress et al.l 
Il992h . The most significant frequencies in the periodogram are 
located at periods ~0.5d and ~840d, as shown in Fig. The 
time series of observations is insufficient to provide a unique 
period and more observations are required in order to better con¬ 
strain any periodic behaviour; however, we verify that the mea¬ 
surements phase well with these periods, and do not appear to 
vary in a coherent way when phased with other periods corre¬ 
sponding to lower peaks in the periodogram. 

Due to the limited temporal sampling of our dataset and the 
precision of our measurements, it is difficult to directly measure 
the expected radial velocity differences from spectra taken on 
the same night. Therefore, it is not possible to verify the plau- 
sibility of the ~0.5 d period solution. However, as shown by 
iNeiner et alJ (|2012|) . rapid radial velocity variations that occur 
over the timescale of a single polarimetric sequence can induce 
detectable signatures in the diagnostic null profiles. These signa¬ 
tures result from residual polarisation signals that were not prop¬ 
erly cancelled during processing because of the radial velocity 


shifts. Therefore, the polarimetric spectra afford us the oppor¬ 
tunity to test for the presence of short-period variations. To do 
so, we compared a synthetic modeQ which takes into account 
the predicted radial velocity variations according to the ~0.5 d 
period solution, to the mean LSD profiles extracted from the 
ESPaDOnS spectrum to test if null signatures should be present. 
The ESPaDOnS spectrum was obtained with the longest expo¬ 
sure time, and therefore should show the largest effect due to 
velocity shifts. Because of the strong polarisation signal, our re¬ 
sults show that when radial velocity shifts are added to each in¬ 
dividual sub-exposure (in agreement with the expected variation 
suggested by short-period variations), the resulting null profile 
should also show a strong signature, which is easily detectable 
in the mean LSD profile. While this test cannot definitively rule 
out the possibility of short-term radial velocity variations, it does 
provide compelling evidence against it. This result suggests that 
the ~0.5 d period is probably an alias of the correct period. We 
tested the plausibility of this hypothesis by subtracting a sinu¬ 
soidal fit to the RV variations corresponding to the ~840 d period 
plus its first harmonic, and recomputed the periodogram on the 
residuals. The resulting periodogram no longer shows any sig¬ 
nificant power about 0.5 d, demonstrating it to be an alias of the 
-840 d period. 

If the radial velocities do vary with the -840 d period, it is 
our conclusion that these variations are likely the result of bina- 
rity. If these variations were due to shifts in the centre-of-gravity 
of the line profile due to the inhomogeneous surface distribution 
(e.g. spots) that is common among Ap/Bp stars, then we would 
expect much smaller radial velocity shifts, which reflect the line 
distortion, with a maximum velocity range of the order of the 
line width. As well, we would expect periodicity consistent with 
the rotational period or one of its harmonics. 

If we adopt this period as the orbital period then we ob¬ 
tain the orbital solution given in Table |7] using ver sion 1.0.2 of 
the pr ogram losp (Liege Orbital Solution Package; [Rauw et all 
l2000h . Eigure |7] displays the RV orbital curve for this long pe¬ 
riod. Of particular interest, this solution suggests a high eccen¬ 
tricity -0.4 and a high mass-function f(m) - m? sin^ il(M + 
mf' - 0.4 (where M is the mass of the observed star, m is the 
mass of the unseen companion and i is the orbital inclination). 
Using the Icr limits of our estimated mass of HD 94660, this 
mass-function implies a lower mass limit for the companion of 
>2 Mq. However, after co-aligning our spectra we find no evi¬ 
dence to suggest the presence of any additional spectral features 
that are not associated with HD 94660, which should be visible 
if HD 94660 hosted a -2 Mq MS companion. Given the rela¬ 
tively high mass-function permitted by our preliminary orbital 
solution, the most likely candidate is a high-mass neutron star 
or black hole; however, further data is required to constrain the 
period and verify the orbital solution before any definitive con¬ 
clusions can be made. Eurthermore, our solution does not rule 
out the possibility of a hierarchical system, where the compan¬ 
ion is a combination of several objects with a total combined 
mass of >2 Mq (such as two white dwarfs). 


' This model is constructed by producing a sequence of individ¬ 
ual synthetic sub-exposures with different radial velocities and treat¬ 
ing the m in the same fashio n as the observations using the double-ratio 
method dPonati et allll997f) . 
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Period (d) 


Fig. 6. Lotnb-Scargle periodogram from the radial velocity vari¬ 
ations of HD 94660. 



Phase 


Fig. 7. Radial velocity orbital solution for the long 840 d pe¬ 
riod (dashed red line). Shown are the RV measurements for the 
HARPSpol (red squares), UVES (blue triangles) and ESPaDOnS 
(black dot) spectra. 

8. Discussion and Conciusions 

HD 94660 is an Ap star commonly used as a magnetic standard 
for polarimetric observations in the southern hemisphere. It has 
an effective temperature Teff = 11 300 K with logL/Lo = 2.02 
and mass M/M© = 3.0. The rotation period is approximately 
2800 d and v sin i is less than about 2 km s '. The surface mag¬ 
netic field strength is of order 6 kG globally. 

The aim of this project is to establish a preliminary mag¬ 
netic field model of HD 94660 to use to estimate the atmospheric 
abundances of several chemical elements. The magnetic held 
model adopted is a simple, low-order axisymmetric multipole 
expansion whose parameters are established by htting the ob¬ 
served periodic variations in {B^} and (B) to computed models. 
This model is produced in the framework of the oblique rotator 
model and reasonably reproduces the observed variations in {B~) 
and (B) with rotational phase (Pig.[Tll. The model is only a coarse 
approximation to the true held geometry of HD 94660, but is 


Table 6. Log of radial velocity variations of HD 94660. 
Recorded are the Date and JD of the observations and RVs. 
The uncertainty in the measurements of RV are less than about 
+ 1 kms^'. 


Date 

(DD-MM-YYYY) 

JD 

(2450000+) 

RV 

(km s“') 

01-05-2001 

2031.464 

16.8 

01-08-2001 

2038.441 

16.9 

03-12-2005 

3707.841 

18.1 

09-01-2006 

3745.167 

16.6 

24-05-2009 

4975.546 

19.4 

25-05-2009 

4976.536 

19.8 

31-05-2009 

4982.605 

21.8 

05-01-2010 

5201.833 

30.2 

19-05-2011 

5701.450 

0.4 

20-05-2011 

5702.449 

0.4 

01-04-2012 

6018.559 

32.2 

28-04-2014 

6775.610 

35.6 


Table 7. Preliminary orbital parameters for HD 94660. Included 
are the orbital period (P), the time of periastron (To; 2450000-I-), 
the centre-of-mass velocity (y), the eccentricity (e), the velocity 
semi-amplitude (K), the angle of the line of nodes (a>), the pro¬ 
jected semi-major axis (a sin /), the mass-function (/(m)) and the 
standard deviation ht (cr). 


Parameter 

Best-fit 

P(d) 

840* 

To 

1611 ±3 

3 

1 

18±0.1 

e 

0.38 ± 0.03 

K (kms“‘) 

17.7 + 0.3 

cj (deg) 

271 ±2 

a sin i (Rq) 

272 ±6 

f{m) (Me) 

0.39 ± 0.02 

cr 

0.41 


* - Hxed 


able to reproduce the observed Zeeman splitting and Stokes V 
signatures with rotational phase reasonably well (Pig. [3]l. This 
simple magnetic model is adequate to make it possible to de¬ 
termine a hrst approximation of the atmospheric abundances of 
this star. The actual parameters of this model are discussed in 
Sect. 4.3 and presented in Table 1. 

We have used a dozen high-dispersion I spectra, well dis¬ 
tributed in phase over the rotation cycle of the star, for a prelimi¬ 
nary investigation of the surface chemistry and a characterisation 
of how the derived abundances may vary over the stellar surface. 
Prom the magnetic field model, the fact that i plus fi is small 

O 

(< 50 ) indicates that we are mainly observing one magnetic 
hemisphere, with limited information about the opposite mag¬ 
netic hemisphere. Although more than half of the stellar surface 
is seen, our investigations indicate very little abundance vari¬ 
ations with rotational phase. Por all elements studied, a single 
abundance hts well all available spectra and therefore a model 
with a uniform abundance distribution over the stellar surface is 
adopted. 

As is expected for magnetic Ap stars, most elements stud¬ 
ied have non-solar abundances. The abundances of O, Mg and 
Ca are all slightly below solar abundance ratios. Only an up¬ 
per limit for He is possible which clearly classifies HD 94660 
as He-weak, with a value at least 10 times less than the solar 
abundance. All other elements studied are more abundant than 
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in the Sun. Most drastically, the rare-earth elements La, Ce, Pr 
and Nd are all between about 10^ to 10“* times more abundant 
than in the Sun. The abundances of Si, Sr and the Fe-peak ele¬ 
ments Ti, Cr, Mn, Fe, Co and Ni are also larger than the solar 
abundance ratios (by factors of order 10^ or less). Although no 
significant variations with co-latitude are observed in the stellar 
spectra, within a single spectrum the Fe-peak elements, calcium 
and silicon exhibit strong evidence of vertical stratification. This 
is most notable for Fe, where abundances derived from weak and 
strong lines of Fe ii (as well as weaker neutral Fe lines) differ by 
factors of order 10. Furthermore, the discre pant abundances de¬ 
rived fr om lines of Si ii and Si iii reported bv iBailev & Landstreej 
(l2013ah are also present in HD 94660. A more detailed analysis 
of abundance stratihcation in this star is clearly warranted. The 
discovery of Co in HD 94660 was surprising and further detailed 
studies of other hot magnetic Ap stars are recommended to as¬ 
certain the overabundance of this element in the atmospheres of 
these stars. 


iLandstreet et al.l (|2014|) highlight the fact that HD 94660 is 
a star where there is poor agreement between the mea¬ 
surements made from instruments with lower and higher res¬ 
olutions. This is evident in Fig. [1] where the held strengths 
extracted from the HARPSpol and ESPaDOnS spectra using 
the entire metallic spectrum disagree with the FORSl mea¬ 
surements. This phenom enon is also presen t in other stars 
such as HD 318107 (see iBailev etall 1201 ih. NGC 216 9-12, 
NGC 2244-334 and HD 149277 (see Bagnulo et al.Ll2006ll . The 


disagreement between held measurements, using LSD, of dif- 
ferent elements found in HD 94660 is also not uncommon (e.g. 
iManfroid & Mathvsl(l2000t) : lBailev et aT](l201lh for HD 318107 
and iBailev et al.l (I 2 OI 2 I 1 for HD 133880). These types of dis¬ 
crepancies are generally considered an indication of the inho¬ 
mogeneous held distributions and large horizontal abundance 
variations (“spots”) on the stellar surface. This hypothesis is 
supported by more detailed maps of Ap stars using magnetic 
Doppler imaging (MDI) in which clear, complex held distribu¬ 
tions and an omalous abundance spot s are observed on the stellar 
surface (e.g. lKochukhov et al.Ll2004l) . 

The RV variations measured in HD 94660 arise from orbital 
motion. Because of the long period currently favoured by the 
measurements (~840 d), this would suggest a massive compact 
companion such as a high-mass neutron star or black hole, or 
possibly a hierarchical system, where the companion is a combi¬ 
nation of several objects with a total combined mass of >2 Mq. 
It is rar e for an A-type star to host a massive compact compan- 


_ ffPe _ 

ion (see lKaner et al.l 120061) and further monitoring is warranted 
to better constrain the properties of the companion. This sys¬ 
tem could help to establish the role that binarity may play in 
the origin of magnetism i n stars with radiative envelopes (e.g. 
iGrunhut & Alecianll20l4 and references therein). 

HD 94660 is a star that warrants further investigation and 
highlights the need to study more sharp-lined magnetic stars in 
detail. Such studies are crucial to further understand the inter¬ 
play between the magnetic held and the formation of vertical 
stratihcation in the atmospheres of Ap stars. They also provide 
important laboratories to test our multi-line techniques for mea¬ 
suring magnetic helds, such as LSD, by allowing measurement 
of magnetic held strengths from individual lines, a task that is 
not possible for stars that are fast rotators. At present, it is un¬ 
clear what the discrepant held measurements for different lines 
of the same element are telling us about the chemical or mag¬ 
netic structure of Ap stars, and therefore further analysis is re¬ 
quired. Long term monitoring is also clearly indicated to hrmly 
establish the nature of the RV variations observed in HD 94660. 
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